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Aerodynamics  Research  Report  No.  271 


FORMULATION  OF  TtiE  COMPLETE  EQUATIONS  OF  BOUNDARY  LAYER 
STABILITY  WITH  MASS  TRANSFER 


by 

John  0.  Powers 


ABSTRACT:  The  stability  theory  of  laminar  boundary  layers  in 
response  to  infinitesimal  disturbances  is  re-examined  for  the  case 
of  a  binary  mixture  with  foreign  gas  injection.  Because  of  the 
inherent  limitations  of  the  asymptotic  stability  calculation  pro¬ 
cedures,  an  approach  was  taken  which  utilizes  the  complete  equation 
system.  Such  an  approach  is  described  herein  and  the  resulting 
system  of  equations  is  presented  in  a  manner  suitable  for  numerical 
evaluation. 
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Formulation  of  the  Complete  Equations  of  Boundary  Layer  Stability 
with  Mass  Transfer 

The  present  analysis  extends  the  methods  of  direct  solution  of  the 
boundary  layer  disturbance  equations  to  account  for  the  effects  of 
foreign  gas  injection  on  the  laminar  boundary  layer  stability.  The 
method  is  based  on  the  work  of  W.  B.  Brown  and  L«  M.  Hack,  with  the 
differences  coming  from  the  addition  of  a  species  continuity  equation, 
diffusive  flux  terms  in  the  energy  equation,  a  modification  of  the 
form  of  the  equation  of  state,  and  a  difference  in  the  dependency  of 
the  transport  properties  on  the  state  variables. 

The  present  report  covers  only  the  formulation  of  the  equation  system 
to  a  point  where  solutions  may  be  sought  by  numerical  methods. 
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f(  TJ) 


On(Tl) 


SYMBOLS 

factor  in  the  general  solution  of  equation  (98)  defined  in 
equation  (129) ,  and  given  by  equations  (130)  through  (137) 

complex  coefficients  defined  by  equation  (60) 

factor  defined  in  equation  (114)  and  given  by  equations  (123) 
through  (128) 

complex  coefficients  which  are  combinations  of  the  a , - ;  given 
by  equations  (103)  through  (113)  J 

specific  heat 

specific  heat  ratio  defined  in  equation  (23) 
mass  concentration 

coefficients  defined  by  equation  (129) 
imaginary  component  of  disturbance  velocity 
real  component  of  disturbance  velocity 
binary  diffusion  coefficient 
molecular  diameter 
strain  tensor 

function  defined  by  equation  (12) 
mode  function  for  u' 
quantity  defined  in  equation  (34) 
enthalpy 

quantity  defined  by  equation  (83) 

length  scale  used  in  non-dimensionalizing  the  equation 

syste">  "  /sir 

X 

molecular  weight 

mode  functions  for  the  variables  given  in  page  8, 
n  -  1,  2,  3,  4,  5,  and  6 

Prandtl  number 
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SYMBOLS  (Cont'd) 

p  pressure 

Q  general  designation  for  a  quantity 

R  gas  constant 

r(rj)  mode  function  for  p* 

Re  Reynolds  number 

Sc  Schmidt  number 

T  temperature 

u  velocity  parallel  to  plate 

u.(or  u.)  general  designation  for  velocity  in  the  itfa  (or  jth) 
J  direction 


V  set  of  dependent  variables  defined  in  page  19, 

n  n  -  1,  2,  3,  and  4 

v  velocity -normal  to  plate 

w(ti)  dimensionless  velocity  ratio,  defined  in  page  8 

Cartesian  coordinate  in  the  i  direction 
y  coordinate  normal  to  the  plate 

Z  set  of  dependent  variables  defined  in  equation  (59), 

n  -  1,  2,  3,  4,  5,  6,  7,  and  8 


a 


n 


8(t}) 

X 


wave  number  of  the  disturbance 
quantity  defined  in  equation  (34) 
ratio  of  specific  heats 
Kroenicker  delta 

Blasius*  similarity  variable,  ^  /Tte^ 

quantities  defined  in  equations  (55)  through  (58) , 
n  -  1,  2,  3,  and  4 

mode  function  for  h 1 

viscosity  coefficients 
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SYMBOLS  (Cont'd) 

X1  exponents  in  the  particular  solutions  (Eq.  (114))  given 

by  equations  (119)  through  (122) 

p  coefficient  of  bulk  viscosity 

v  coefficient  of  kinematic  viscosity 

5 (q)  mode  function  for  C^' 

n(?l)  mode  function  for  p / 

p  density 

t  time  scale  used  in  non-dlmensionalizing  equation  system, 

t 
u 

stress  tensor 

qp(H>  mode  function  for  v' 

Subscripts 

a  refers  to  mean  flow  gas 

b  refers  to  injected  gas 

<5  refers  to  boundary  layer  thickness 
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INTRODUCTION 


The  effect  of  foreign  gas  injection  on  the  stability  of  the 
laminar  boundary  layer  was  first  analyzed  by  Shen,  reference  (1) , 
in  1957  using  asymptotic  procedures.  He  developed  the  "inviscid 
solutions"  in  power  series  of  the  square  of  the  wave  numbers,  a3 , 
and  then  demonstrated  how  the  "viscous  solutions"  could  be  obtained 
by  an  asymptotic  expansion  in  the  inverse  square  root  of  the  product 

of  the  wave  number  and  the  Reynolds  number,  (a  Re^)”^,  By  further 
manipulation  he  was  able  to  show  that  the  influence  of  foreign  gas 
injection  wa3  evident  primarily  through  the  mean  boundary  layer 
profiles  and  through  a  simple  correction  factor  in  the  well-known 
Dunn  and  Lin  secular  equation,  reference  (2) .  In  additicn,  Shen 
developed  an  "inviscid  criterion"  for  the  stability  of  injection 
profiles  and  by  applying  this  criterion  he  was  able  to  demonstrate 
that  the  injection  of  a  heavy  molecular  weight  gas  as  a  coolant 
might  lead  to  boundary  layer  with  improved  stability  characteristics. 

The  possibility  of  improving  the  stability  characteristics 
of  boundary  layers  by  injection  of  heavy  molecular  weight  gases 
stimulated  the  research  of  Powers,  Heiche,  and  Shen,  reference  (3), 
who  made  a  qualitative  investigation  of  this  phenomenon.  They 
reformulated  the  asymptotic  procedures  of  reference  (1)  to  facili¬ 
tate  numerical  solution  of  the  stability  equations  and  made  a 
systematic  investigation  of  the  effects  of  varying  the  molecular 
weight  and  diameter  of  the  injected  gas.  It  was  shown  that,  in 
terms  of  minimum  critical  Reynolds  number,  the  injection  of  a  small- 
diameter  i ight-weight  gas  could  decrease  the  stability  by  as  much 
as  an  order  of  magnitude.  In  contrast,  the  injection  of  a  large- 
diameter  heavy  gas  could  actually  improve  the  stability  of  the 
boundary  layer.  Those  results  were  for  zero  Mach  number  thermal 
boundary  layers;  however,  the  findings  were  found  to  apply  up  to 
low  supersonic  Mach  numbers  in  a  later  investigation,  reference  (4) . 

Attempts  to  extend  the  results  to  Mach  numbers  higher  than  about 
1.3  by  the  asymptotic  approach  were  not  possible  because  of  an 
inherent  limitation  in  the  procedure.  This  limitation  is  believed 
to  be  associated  with  the  apparent  singularity  in  the  "inviscid" 
equation.  The  singularity  occurs  when  the  velocity  of  the  distur¬ 
bance  relative  to  the  wall  becomes  supersonic  and  as  a  result  one 
of  the  coefficients  of  the  differential  equations  changes  sign.  This, 
in  turn,  violates  one  of  the  conditions  of  the  Sturm-Licuville, 
theorem  which  guarantees  the  existence  of  eigensolutions  only  under 
certain  conditions.  It  has  been  possible  to  obtain  solutions  at 
high  Mach  numbers  by  changing  dependent  variables  in  a  manner  suggested 
by  Lees  and  Lin,  reference  (5) .  Such  solutions  were  obtained  by 
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Reshotkc,  reference  (6),  but  they  were  not  unique  eigensolutions. 

L.  Mack,  reference  (7),  demonstrated  in  fact  that  the  multiplicity 
of  solutions  increased  as  the  Mach  number  increased. 

The  use  of  asymptotic  procedures  at  high  Mach  numbers  would 
require  the  inclusion  of  many  terms  which  had  previously  been 
neglected  and  as  a  result  the  numerical  program  would  become  very 
complex.  Under  these  circumstances  it  was  considered  desirable 
to  investigate  the  high  Mach  number  effects  of  mass  transfer  on  the 
boundary  layer  stability  by  developing  a  method  of  directly  solving 
the  complete  linearized  disturbance  equations  by  numerical  methods. 
The  feasibility  of  the  direct  solution  method  has  already  been 
established  by  the  excellent  pioneering  work  of  Drs.  W.  Byron  Brown, 
reference  (8),  and  L.  Mack,  reference  (9).  Both  of  these  investi¬ 
gations  have  been  well  documented  and  may  be  considered  as  the 
foundational  work  upon  which  the  present  development  is  based.  It 
is,  therefor^  the  objective  of  the  present  study  to  extend  these 
previous  methods  to  accommodate  the  effects  of  foreign  gas  injection 
on  the  laminar  boundary  layer  stability.  The  present  investigation 
is  a  formulation  of  the  equation  system,  and  methods  of  obtaining 
numerical  solutions  are  indicated.  Numerical  solutions  are  not 
included  in  the  present  report  since  it  was  considered  desirable  to 
expedite  the  dissemination  of  the  formulation  before  complete  solu¬ 
tions  were  obtained. 


ANALYSIS 


General  Equations 

The  general  equations  of  motion  for  a  binary  gas  system  in 
Cartesian  tensor  notation  are  used  as  the  basic  equations  for  the 
present  analysis.  These  equations  are  the  equations  of  motion, 
global  continuity,  species  continuity,  and  energy.  In  the  dimen¬ 
sional  form  they  are  the  following: 


Motion: 


Global  Continuity: 


Species  Continuity: 


>f 

1  Jt  4 


(i) 


(2) 


(3) 
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^  i  2^-J-  at 


,  2>p  (-=■ 

+  U\  —  +v- 

J  ^  *, 


£ 

+  St; 

V 


(4) 


In  these  equations  all  quantities  take  their  conventional  definitions 
as  defined  in  the  list  of  symbols.  The  terms  which  result  from  the 
inclusions  of  mass  injection  are  related  to  the  quantities:  Cj,,  the 
mass  concentration  of  the  injection  gas;  £aj,,  the  binary  diffusion 
coefficient  for  the  diffusion  of  gas  "b"  into  gas  "a";  Sc,  (p/p£>) , 
the  Schmidt  number;  and,  ha  or  h^,  the  enthalpies  of  the  components 
species.  The  stress  and  rate-of-strain  tensors  take  their  usual  fora: 


where  the  viscosity  coefficient,  X,  is  equal  to  three  halves  of  the 
bulk  viscosity  coefficient. 

In  addition  to  equations  (1)  through  (4)  we  must  add  to  the 
basic  set  a  form  of  the  equation  of  state  which  is  compatible  with 
our  application  to  the  binary  system.  This  is  accomplished  by  using 
the  conventional  relation: 


fXT 


(7) 


where  for  the  binary  mixture  the  assumption  is  made  that  the  gas 
constant  and  the  mean  specific  heat  are  functions  only  of  the  species 
concentrations.  This  leads  to  the  relations: 


(8) 
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(9> 


and 


7=  */ej> 


(10) 


By  using  equations  (8),  (9),  and  (10)  in  equation  (7)  and  forming 
the  logarithmic  derivative  the  equation  of  state  takes  the  form: 


where  - 


Cpa.  *-C6  CCfy-C/ >m) 


(11) 


(12) 


In  equation  (11)  the  prime  is  used  to  indicate  the  fluctuation  or 
disturbance  part  of  a  quantity  which,  when  added  to  the  time  inde¬ 
pendent  basic  flow  (indicated  as  barred) ,  yields  the  instantaneous 
value  of  the  quantity.  This  is  written  as: 


%  *  4(*;)  *-  Q'Ctyt) 


(13) 


(Note:  Later  the  primes  will  be  used  to  indicate  derivatives  with 
respect  to  the  independent  variable.  This  notational  change  is  for 
convenience  and  wixl  be  noted  when  it  occurs.) 


Introducing  the  instantaneous  form  of  the  variables  equation  (13) , 
into  equations  (1)  through  (4) ,  subtracting  the  mean  flow  equations, 
and  neglecting  quadratic  terms  in  the  disturbance  quantities  results 
in  the  general  form  of  the  disturbance  equations  as  follows: 


Motion: 


r* 


at 


9x ; 
v 


(14) 
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Global  Continuity:  ^  f  ^ j  **"  f  ^0—  O 


(15) 


Species  Continuity 


:  fhl 


_  2. 


2<  -  2)(V  ,  2C,  7  . 
3 1  1  SXj  J  2>x;  j  1 

f<2) 

2x,-  ^  at. 


» - 

+  PU.-  — 

r  <  >xj 


(16) 


Energy:  p'iL  —  +  P  T  —  4.  U-  —  +  U-C-  — 

1  J  2>Xj  I  (St  J  SK-j  tl)  2>IC 


2*' 

=  ir 

2)  t 


+  9 


CL  f  u/  |£  a'  +<„  fi*.  +- 

9  f  T  *l'h(^'2£L  £]>/> 

32p2}  ax/^UJ  2a,-  r  3/  /  j’c  C' 


‘  j  T 


(17) 


+gt6-yo* 

)l 


[i  C"  *. >]  Xyf)  *  f  i 


r>„  / 


Two  Dimensional  Parallel  Flow 


The  only  simplifying  assumptions  are  now  made;  namely,  that  the 
mean  flow  is  parallel  and  that  the  disturbances  and  the  mean  flow  are 
two  dimensional.  The  first  of  these  assumptions  implies  that  v  «  u 
and  that  derivatives  of  mean  flow  quantities  with  respect  to  x  may 

be  neglected,  (  i.e„,  «  -^p  .  The  second  assumption  eliminates 

the  coordinate  dependence  from  the  problem.  These  assumptions  are 
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conventional  in  boundary  layer  stability  and  no  attempt  is  made  here 
to  evaluate  their  significance  since  this  has  been  the  subject  of 
previous  investigations.  (See,  for  example,  reference  (2).)  As  a 
result  of  these  assumptions,  equations  (14)  through  (17)  can  be 
written  as  follows: 


x-momentum : 


*3 


a£u' 

u  7>jC 


b 


Cbu)  l 


y 


3x2;  >y 


9£  , 

2<i  L 


•dr's. 
2*  J 


+  * 
7  0 


(18) 


Global  Continuity: 


A  , 

A  C)t/'  - 


(20) 
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Species  Continuity;  f  ^ 

r  L 


+  iZd<±\  ..,}&>? 


ri  +  y'f* 

9**0  ~  x-  ^ 

V  hfrht'iy  (-£)  °£-Y£  °  rsi  /  I!u 


Energy ; 

/r^'  Or 


5 ^  +  *2*  SL  (yO  f 

T  7  8 /  L  sJ 


?[***£"•% 7- 

o/„  ,l*  tjJ 


L  ,  /  /  cJCl  x.  0/-  „  J 

^  (^^9f  (%'+¥)  +  *£ >C%  ^ 
n/  7  V  V  <22> 


2£L  2l'  1  ,  jt  r  /  - 
?/  /  Sj  Oj  76  ^ 

'  7  1  ?  &  7—7: - ~T  "1 


7‘  ~  y/'  f 

,  -7  24 '  >-  ,  r  Q  <r  x 

x  r  a7'  f '■*%<[*(%*%, 

'  *'  £*  $  'S/Wg',  *',  _ ; 


r-S'. 


l  _  1  &>/  <7 

*  "A  ^  <%  z—z'^'K,'  nr~ 

*  sr  >rl»**S])*f(%+f)]h 


)  ^ 


)i  *l  $ 


1  °} 


Cf(  %>  2J~ 

>/  SJ 
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where : 


(23) 


The  Mode  Function  Equations 

It  is  next  desirable  to  express  the  equations  (11)  and  (18) 
through  (22)  in  dimensionless  form.  This  is  accomplished  by  intro¬ 
ducing  the  length  scale,  i  m  x//Rex  ,  the  time  scale,  r  -  i/u^  ,  and 
by  scaling  the  other  variables  with  respect  to  their  values  at  the 
edge  of  the  boundary  layer.  The  concentration  of  the  injectant  is 
an  exception  to  this  scaling  since  it  is  already  dimensionless.  We 
note  that  scaling  results  in  the  use  of  the  Blasius  variable, 

r\  -  —  /Re  ,  to  represent  the  y  coordinate.  It  is  further  desirable 

A  A 

to  represent  all  of  the  fluctuations  as  harmonic  functions  whose 
amplitude  is  determined  by  a  "mode?*  function  of  q.  Using  these  con¬ 
siderations,  the  mean  flow  and  fluctuation  quantities  are  expressed 
as  in  the  following  table: 


Mean  Flow  Variable 


TABLE  I 

Fluctuation  Quantity 


u,  -  U)  Cj) 

1/  ~  Uoc  V' Oj') 

^  '=u0o( 

f-  fco  fcp 

P1-  £  Hy)  £ '*  Cx~c^ 

P  --  A  P 

p'=  ^  rrCij)  a  Cx'Ct) 

Pm  POl) 

j 

II 

« 

eh'=  PCyj)&^Cx-ct) 

A  =  A()j) 

>  '  = ^  Qttflj)  6  (*'ct) 
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w 

7° 


; 


TABLE  I  Cont'fi. 

.A 


c<*  Ck  -  ctj 


-t 


%  A 


cc<Cx-ciJ 


A  ’^r)  --  t^io 


c-^fe-  &>'  = 


r^/=  c%L\A<i^ce> 


Using  the  definitions  given  in  Table  I,  introducing  the  Reynolds 
xj  X  „ 

number,  R  -  — ~  -  /Re  ,  and  using  the  global  continuity  equation 

^  «D  ^ 

in  the  energy  equation  gives  the  desired  form  of  the  mode  function 
equation.  (Note:  For  the  following  pages  the  primes  will  be  used 
to  indicate  derivatives  with  respect  to  the  variable  Tj.) 


x-momentum 


:  f  ^  "  A  if  '-if)] 


SCS/'A+aZ  fay*  %'*+>'(/«*)’)} 


(24) 


y-raomer.tum;  l  (lO-C) 


-f7J 


Jl  ,r/  V 


7  if  1  f"f')  *it  {<*/■>*//, | 


(25) 


+-  0 


9 


''•s- 


Global  Continuity: 


vt  at  mn 


n  wjj  ±  n. 


/?/?  1  0«7 

uu-i.o  t 


!  (v-c)  r  y-  P  c  f  +<f')  -O 


(26) 


Species  Continuity: 


f  +- <Vf  }  =  (27) 

Energy :  j f>  [i  Cu>-C)  rf'J  ~  ~  jd/C^-OTT  -  !&<-)  [  f 

+ ft i+  (28> 

-f/*  y.  -/(f-4)  ¥fy[(i-^  ietl 


Equation  of  State : 


±.  £  ^  Jl_ 

7F=  ^  a.  F(Ct) 


(29) 


Before  equations  (24)  through  (29)  can  be  developed  further  it  is 
necessary  to  establish  the  form  to  be  used  for  the  thermodynamic  and 
transport  properties  of  the  mixture. 

Thermodynamic  and  Transport  Properties 

One  of  the  major  differences  of  the  present  development  and  the 
previous  developments  of  Brown  and  Mack,  references  (8)  and  (9),  is 
the  manner  in  which  the  thermodynamic  and  transport  properties  enter 
the  equation  system.  In  the  previous  works,  these  properties  are 
considered  to  be  these  of  a  single  species  and  hence  are  only  functions 
of  temperature  £c  given  pressure.  In  the  present  investigation, 
since  we  are  dealing  with  the  injection  of  a  foreign  gas  into  a 
boundary  layer,  these  properties  are  necessarily  function*  of  both 
temperature  and  concentration.  Since  enthalpy  is  used  herein  in  place 
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of  temperature,  the  fluctuations  of  these  properties  and  their  r\ 
derivatives  will  be  related  respectively  to  the  enthalpy  and  concen¬ 
tration  fluctuations  and  to  the  enthalpy  and  concentration  profiles. 
In  the  mode  function  equations,  the  properties  appear  both  as  r\ 

derivatives,  and  as  their  fluctuation  mode  functions,  m^T}). 

Because  of  the  dual  dependence  these  expressions  are  written  as: 


^<p/  _  dl/ 


vf/J 


and: 


(30) 


(31) 


where  subscripts  h  and  indicate  derivatives  with  respect  to  those 
quantities. 

Further  development  of  these  expressions  is  dependent  on  the 
choice  of  thermodynamic  and  transport  property  relations.  To 
achieve  a  reduction  in  complexity,  the  forms  used  by  Korobkin, 
reference  (10),  are  used.  These  mixture  properties  only  appear  as 
dimensionless  ratios,  with  the  exception  of  the  P.andtl  number  and 
the  ratio  of  specific  heats  at  the  edge  of  the  boundary  layer.  These 
two  then  become  the  pure  air  values.  The  dimensionless  ratios  of 
the  mixture  properties  are  developed  in  terms  of  the  temperature 
ratio,  T,  concentration,  C5,  molecular  weight,  ma  or  mb,  and  molecular 
diameter,  da  or  db,  to  be  the  following: 


i 

lar*/  c  d6J 

/-A  Ct) 


“>  (-  /-Ci  ) 


(32) 


l-f 


I  "A  C%) 


"A 


(33) 


where 


=  ( 


a  r  ava  t2- 


f  ^>4  1 

lint 


(34) 


il 
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9-  r*- 


(35) 


In  addition  we  write  for  convenience: 


#6  -  /f«- 


(36) 


Using  these  expressions  in  the  proper  combinations,  the  required 
relations  in  equations  (30)  and  (31)  become: 


ll 

(37) 

(38) 

^/6  " 

-<&■  +JI 

oZ  '"'w 

(39) 

II 

£ 

(40) 

Cto,c 

j-C- 

(41) 

£ 

o 

II 

-Cf- 

(42) 

- 

°?e/>3 

(43) 

u 

j*  C^k ) 

(44) 
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The  equation  system,  equations  (24)  through  (29) ,  is  reduced  to 
normal  form  by  the  introduction  of  a  new  set  of  variables  defined  as 
follows: 

2/*^;  ***/■';  *3-f;  ;  Zr--0  <59> 

--  6  ;  2f  =  ?  ;  4  =  / / 

In  this  form  we  have  a  system  of  eight  first  order  differential 
equations. 

* 

2v  -  Z  7^'  ( l  =  /,...  g)  (60) 


The  coefficients,  a^j,  are  made  up  of  combinations  of  the  mean 
velocity,  enthalpy,  density,  and  concentration  profiles,  the  thermo¬ 
dynamic  and  transport  properties,  and  the  four  basic  parameters  of 
the  problem;  namely,  the  wave  number,  a,  the  Reynolds  number,  R,  the 
real  component  of  the  wave  velocity,  cr,  and  the  imaginary  component 
of  the  wave  velocity,  ci.  Of  the  64  possible  complex  coefficients, 
j ,  only  36  are  nonzero  and  they  are  the  following: 
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- 

/ 

(61) 

- 

cf*(jeC»-z) 

~^z - ^°< 

(62) 

aa  = 

--i-  _  /O'!,  \  p ' 

04  ^  C  ^  ^  /  c<» 

(63) 

a„  — 

oC'  z^7  O^  -.  ^  o  ^  1 

^  /'A-  ^  <5  /  ^/""Z^j 

(64) 

a»  - 

'•"S  a 

3  O  c-)(/*-J^  ) 

(65; 

J7  0*-^)-^ f»-sfcu  ,  f)j 

(66) 

tZ-ll.  = 

J.  v' 

•z  ^ 

(67) 

- 

^  "J  " &-£)('*•  (Cp  +&--%') 

(68) 

^Ig  - 

k/'(  ^ 

«<  ^4*  ' 

(69) 

&}l  =  - 

o 

-  C 

(70) 

= 

-r/f 

(71) 

-  ~/°  fcoM#  (l4/-c) 


15 


NOLTR  68-187 


=  /-£.  Cw-c) 


(73) 


3  3}  =  C  S)Ck/'C) 


a 


y/ 


a 


is.  - 


<**  fx.  ,  „  A 

-9l. 


P's 


r 


(74) 


(75) 


(76) 


ci</3 


—  sL 


\%  0*+j;)U  (’f*)  -  ^ >*Z2  J  *  a*\ 


l 


(77) 


a 


a  - 


-  i-  iJiJS  (J^u  C„.c j./c^,'j<"> 


^itr 


3 

a 


*T  /  /  J5y:  — \  I  0  /  <fod  ,  ~{,  ~L  k/ '  1 

s  x  r*  i  <^«J  /  z] 


a 


i(o  ' 


J*jL  /  ^  / 

<&*  ^tSy'  * ) 


(79) 


(80) 


^£ 

j/ 


ffcpjpty-t)- (“<>[(  ?*■**)(  Cfti) 


Q.js  = 


(82) 
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where : 


(i=  Jr*  +  ft  LK  c**j;  x*-*# 

/ 


&C2-Z 


~  1 


-±\  c MJCsCp"' 

&Q  /  / 


^2^3  " 


£ 


6^ 


pif** [*■-<>* %'- 

-  «?  *  C  4. -/)  M*  <*/*  *  '  j 

.  „  P, 

r.  -a:  7L  4#  (&->)  s  CW'C.) 


4-  £L1 
IT 


(83) 

(84) 

(85) 

(86) 


(87) 


a*** 


&&,  “ 


£4,  ^  Kftr?/,, 

-  a-f^  f  *  u  >]  j  g  -#4'v 

^4“A'V«'  f 

,  $  /■■  A  «U«4 't  <  w-fot]  -  *•  u-,] 
(A-  f  )-A«"'  &,)'&%' 
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(EQ.  (90 ^  continued) 

r> 


J-^-3  +  ~£'  +<~6  tfyC&C/O- 

A  CfJl3  ")  +  (yiy*3  c/Uij  /-  ~  -C^Cr- 


) 


f  ] *  cp  c*  "£. [fyf+ji]- 3*$  Vjij  ]  *■ 


/  ^ 


//  ^  ^ 
/•  — 


at 


r^fe  V  f-V- 


Cp  (h  C/> *-  3Jl3  *  o< */)  +J\^  /-  £ cp~[3£ 

P)-C’/] 
a“  ’J  {ef  eH[Jl3 


a 


+  f?  [3~i'-s-ep  ep  yC] 

-  i 


(91) 


n 

&$s  -  *zf>  #»«£€* 


Ji 


(92) 


(93) 


a* =  «*  fe*'fr 


rtg<o 


&z?  = 


^  -%L 

c  Pr*  ^ £ f[u<)  ~ r/i-  C  r{r)0-b"  +  ^ ) 

ab/ZCty  Cf^+GrCf  fV 


(94) 


(95) 


(96) 
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J  *  / 


Method  of  Solution 


(97) 


While  the  present  report  does  not  go  beyond  a  formulation  of  the 
equation  system  it  is  possible  to  draw  upon  the  work  of  Brown  and 
Mack,  references  (8)  and  (9),  and  to  develop  a  suggested  method  of 
solving  the  equation  system,  (60)  .  In  both  of  these  references  the 
authors  constructed  linearly  independent  solutions  which  were  inte¬ 
grated  across  the  boundary  layer  and  were  then  combined  to  yield 
a  general  solution  satisfying  the  boundary  conditions.  The  two 
methods  differed  in  the  direction  of  integration  and  in  the  search 
procedure  used  to  vary  the  initial  guesses  until  a  proper  set  of 
eigenvalues  was  obtained.  The  present  suggested  procedure  follows 
more  closely  the  method  of  Mack,  reference  (9). 

Initially,  we  look  for  four  linearly  independent  solutions 
which  are  applicable  at  q  >  q^.  These  solutions  are  found  by 
solving  the  equation  system: 

°°  * 

2?7r  2  Q .  2?-  Ct-'/'-?)  (98) 

;=/  V  J 


where  the  superscript  asterisk  on  the  a^  coefficients  means  that 
they  take  on  their  constant  value  at  q  >  To  facilitate  the 

derivation  of  these  solutions  the  dependent  variables  are  changed 
by  the  substitutions:  -  Z;  "  Z4’  V3  Z5*  aud  V4  "  Z7*  Th:*-S 

changes  the  eight  first  order  equations  into  four  second  order  equa¬ 
tions  which  are: 


H "  = 

^12,  2.  +*  1 

fc,3  V3  4-  b,^  V, 

(99) 

v;, 

Vi  4- 

4-  bn  \ 1  yj 

(100) 

- 

V  2.  + 

b?3  Vj  4-  \J 

(101) 

w 

M »* 

\A| 

(102) 

where  the  b^j’s  are  the  following  combinations  of  the  ajj’s: 
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K  =  a*  a*  (103> 

b(1.  =  a,t  ah  <l04) 

bu  =  0,1  air  (105) 

b*,H  *  ^12-  (106) 

bzj.  -  <3  Q  ^  -t  Qm  (Xi^  Ci^b  <3<rH  (i°7) 

bzj  -  (*W  z  &ir  +-&h<*  Qur~  (108) 


i)iq  = 

+-  -t- 

(109) 

bji.  = 

gA  al. 

(110) 

4  - 

Gk  air 

dll) 

^iy  - 

r  *  + 

&n  <%(,  / 

(112) 

(113) 

The  system  of  equations  (99)  through  (102)  has  particular 
solutions  which  may  be  written: 

/•',  ^V-'W  d 

‘  j  =  8 

Putting  the  "j**1"  particular  solution  into  the  equations  yield 
the  following  equations: 

(b„  -  A/)^(aV  b,2  b,3  B3‘JJ+  b,y  ^  -  O 

Aw  l^l  =  0 


(114) 


(115) 


(116) 
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632  B7  +  (.633  -  J/J  S3  +■  ~bh  1  -  o  (117) 
Cjbsfxj  -  Jj’)  ”B4  u  1  .=  o  (us) 


The  characteristic  determinant  yields  four  negative  roots  and 
four  positive  roots.  The  positive  roots  are  discarded  since  they 
would  give  solutions  increasing  exponentially  in  T).  The  negative 
roots  are : 


X,  =  - 

hi'  -  pij 

hz  -  -{i  Cbu  +  ki)  -  [-i,  Cbji-  bsi)  *■  ihj  kit.]  lj'i' 


(119) 

(120) 

(121) 


X 


(122) 


Now  by  selection  of  a  specific  magnitude  for  one  of  the  Bp^’s 
in  each  of  the  four  j  groups  of  fundamental  solutions  the  equa¬ 
tions  (115)  co  (118)  can  be  solved  to  yield  the  following  values: 


For  j 


0) 


B 


(123) 


y  l 
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For  j  -  2  or  3: 


A •x-  -  bn 

Q 


U\ 


Bjg'-  -ij>- 

=  o 


(124) 


For  j  -  4: 


™  }  C  iH  )"*  ^*3  ^32.  J  “  b>  tl 

b  23  +  i>»3  **  kjH  )-  bii  b51]  ]•  [(bn 

«i  ■ 

[C  bn-  bss  )Cbt3-  bn*')  ~bis  b^]  J 


(125) 


%  -  jf^X^)  ~  b  Chr  t*v)j Jcki-bii)(bji-bii) 

—  bis  b-$h  \  (i 


-l 


J  u26) 


B3  r  bn~  by  (bii~hy)j  (bzi-b^)  (127) 


(128) 


Next,  it  is  possible  to  write  a  general  solution  of  equation  (98)  as: 


n 


1 


/  *  =  c— 

*“j  vJ  t 

vP< 


c-  acVj(7~^) 


(t  =  I;  ••  0)  <129> 
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where  we  have  used  the  fact  that  the  four  exponentially  growing 

( i  N 

fundamental  solutions  are  discarded.  The  A  ’s  of  equation  (129) 

/  •  \  * 

are  now  related  to  the  bH'*s  as  follows: 


■s/1’ 

(133) 

A*  - 

BS'J> 

(134) 

A?  * 

1  h  U> 

(135) 

Af'“  - 

T>  ti) 

(136) 

a 

*i  s<"’ 

(137) 

The  four  j  values  now  yield  four  fundamental  solutions  which 
may  be  numerically  integrated  from  ti  =  ri^  to  the  wall,  n  «■  0.  At 
the  wall  the  fundamental  solutions  are  combined  to  give  the  general 
solutions  which  satisfy  the  remaining  boundary  conditions  for  high 
frequency  fluctuations;  namely, 

2(  =  t3  -  =2^  =  0  U38) 
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CONCLUDING  REMARKS 

The  present  analysis  extends  the  methods  of  direct  solution  of 
the  boundary  layer  disturbances  equations  to  account  for  the  effects 
of  foreign  gas  injection  on  the  laminar  boundary  layer  stability. 

The  similarities  between  the  present  method  and  the  methods  of 
Brown  and  Mack,  references  (8)  and  (9) ,  are  obvious  and  do  not 
require  further  comment.  The  major  differences  in  the  present 
netaod  which  result  from  t.he  inclusion  of  mass  transfer  come  from 
the  addition  of  the  species  continuity  equation,  the  diffusive  flux 
terms  in  the  energy  equation,  a  modification  of  the  form  of  the 
equation  of  state,  and  a  difference  in  the  dependence  of  the  trans¬ 
port  properties  on  the  state  variables.  The  present  report  covers 
only  the  formulation  of  the  equation  system,  however,  the  analysis 
has  been  carried  to  a  point  where  solutions  may  be  sought  by  numerical 
methods . 


< 
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